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We investigate temperature effects in a Fermi gas with imbalanced spin populations. From the 
general expression of the thermal gap equation we find, in weak coupling limits an analytical expres- 
sion for the transition temperature Tc as a function of various possibilities of chemical potential and 
mass asymmetries between the two particle species. For a range of asymmetry between certain spe- 
cific values, this equation always has two solutions for Tc which has been interpreted as a reentrant 
phenomena or a pairing induced by temperature effect. We show that the lower Tc is never related 
to a stable solution. The same results are obtained in strong coupling limit. The thermodynamical 
potential is carefully analyzed to avoid the consideration of the unstable solutions. We also obtain 
the tricritical points for the chemical potential and mass imbalanced cases, and beyond these points 
we properly minimize the thermodynamic potential to find the stable and metastable first order 
transition lines. 
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The recent advances in experiments with ultracold 
fermionic atoms have provided the possibility for the un- 
derstanding of superfluidity in several physical situations, 
from high temperature superconductivity to the pairing 
of quarks in the cores of neutron stars. 

When a two fermion species system have the same 
number of spin- up and spin-down particles, its ground 
state is described by the well known Bardeen-Cooper- 
Schrieffer (BCS) theory of superconductivity [1]. In this 
case, if the temperature is below a certain critical tem- 
perature (Tc), fermions with opposite spins interact near 
their common Fermi surface resulting in pair formation, 
even at arbitrarily weak coupling. For temperatures 
above Tc the system is found to be in the normal state 
i.e., the fermions are unpaired. 

The possibility of pairing formation in imbalanced sys- 
tems, where the chemical potentials or the numbers of 
fermionic species are different, was first pointed out more 
than 4 decades ago by Sarma [Ij. Since then, various 
exotic phases have been proposed for the mismatched 
case, such as the Larkin, Ovchinnikov, Fulde and Fer- 
rel (LOFF)-phase [|, the breached pair superfluid phase 
(BP) [1], deformed Fermi surfaces [5|, phase separation 
in real space @, 0], and also the possibility of new co- 
existing phases in the BEG regime 0. Other pairing 
mechanisms beyond BCS, such as P-wave superfluidity, 
have also been investigated . Observation of phase sep- 
aration between a fully paired superfluid core surrounded 
by the unpaired excess atoms, have been reported inde- 
pendently by the Rice [TO, IH and MIT [H, [H groups. 
This phase separation can be viewed in terms of three 
phase transitions of different nature. Differently from 
the standard BCS thermodynamical phase transition we 
mentioned above, in an imbalanced system the phase 
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transitions that may happen are: I. At zero tempera- 
ture (T), increasing the chemical potential or number 
asymmetry the system undergoes a (first-order) quan- 
tum phase transition to the normal state Q; II. Still at 
T = 0, first order phase transitions occur from the nor- 
mal to a phase separation (PS) phase, and from PS to 
a spatially homogeneous (magnetized) superfluid as the 
interaction parameter l/kpa is varied [1, HJ, with a tri- 
critical point sitting in the fully polarized line and 
III. Lowering the temperature, a system with different 
and fixed number particles phase separates into a unpo- 
larized superfluid core surrounded by a polarized normal 
shell, provided the asymmetry does not exceed a critical 
value [T2I , ll^. Thus, it is of crucial importance to ad- 
dress the issue of the temperature effects in imbalanced 
Fermi gases [l5|. Some recent works have investigated 
the relevant problem of imbalanced Fermi gases at finite 
temperature [H, |T3, [3. However, there are only few 
recent theoretical studies exploring, concomitantly, the 
temperature effects and the stability of pop ulation and 
mass imbalanced Fermi gases BE Sill. 

In this work we investigate thermodynamical phase 
transition in a two species homogeneous system. We 
carry out a quantitative study of the lowest energy state 
of an imbalanced Fermi gas on the BCS side of the res- 
onance. In this regime unexpected temperature effects 
manifest, such as the appearance of two solutions for the 
mean- field (MF) transition temperature < Tc,i < Tc,2- 
A gap A would emerge at Tc,i, increase up to a maximum 
value, and then decrease, vanishing at Tc,2- A premature 
interpretation of this nonmonotonic behavior for A and 
its two critical temperatures is that temperature favors 
pairing 0. On the contrary, we find that temperature 
(heat, to be more precise) always acts in the direction of 
disrupting the fermion pairs, mainly in imbalanced sys- 
tems, at least within the model under consideration. We 
show systematically, employing MF theory, both at weak 
and strong coupling, that the lower critical temperature 
is related to an unstable solution of the gap equation. 
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and thus does not correspond to a true lowest energy 
(stable) state. However, we point out that the reentrant 
phenomena could manifest when the BEC-BCS crossover 
theory is considered. In this case, fluctuation contribu- 
tions (necessary for the introduction of the concept of a 
pseudogap) are taken into account [23I. [2^. 

More important than only identifying the unphysical 
solutions of the gap equation, we demonstrate: (i) Ana- 
lytically we derive, in the weak coupling limit, an expres- 
sion for the second order phase transition line from which 
we found the tricritical points {Sjj.c^Tc) and {Snic^Tc). 
(ii) Numerically, we obtain the physical second- and first- 
order phase transition lines separating the superfluid and 
normal phases, both in the weak and strong coupling lim- 
its, (iii) We find the metastable first order curves in the 
phase diagrams T vs. Sfi, and T vs. (5m, ignored in all 
previous analysis on temperature effects in mass imbal- 
anced Fermi gases. 

The Model. To begin with, let us consider a nonrela- 
tivistic dilute (i.e., the particles interact through a short- 
range attractive interaction) cold fermionic system, de- 
scribed by the following Hamiltonian 



= H 



k.a 



(1) 



k,k' 



where a^, are the creation and annihilation operators 

for the a particles (and the same for the b particles) and 

are their dispersion relation, defined by = — Ma, 
,2 

with = 2^ and /i^ being the chemical potential of 
the (non-interacting) a-species, a = a, 6. To reflect an 
attractive (s-wave) interaction between particles a and b 
we take — ^ < 0. 

The Thermodynamic Potential and the Thermal Gap 
Equation. From the Hamiltonian ([T]) one obtains the 
grand potential from which all thermodynamical quan- 
tities of interest can be obtained: 



^=^ + T.H-Ek-THe- 
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Ek ± are the quasiparticle excitations. 



with Ek = ye^^ + A2, and = ^^f^. Deriving the 
equation above with respect to A, we obtain the gap 
equation 



9k), 



(3) 



where fk and are the Fermi distribution functions 
fki9k = l/(e^^^ + 1), with P = 1/T, where we have set 



s^i/A =1.4 
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FIG. 1: The function /(T/Ao) for different chemical potential 
asymmetries. For small asymmetries 6fi < 6fio, there is only 
one solution for Tc. For S/j < 6/j < 6/j.c there are two solutions 
for the transition temperature. When 6/j > dfic there are no 
solutions for Tc from the gap equation. 



the Boltzmann constant equal to one. The critical tem- 
perature Tc is, by definition, the temperature at which 
A = 0. Then Eq. ^ becomes 



k 



(4) 



where = l/T^. 

Analytical Solutions. In the weak coupling limit, it 
is possible to obtain a compact expression for Tc from 
the equation above. After the momentum integration, 
Eq. can be written as [lH 



Tc 



crAo 
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where a = is a dimensionless parameter reflecting 

the mass asymmetry, M = is the reduced mass, 

Ao is the T = BCS gap parameter in the weak coupling 
limit, 3^(x) = ^(^ + f ) + ^(^ - f ), with ^ being the 
digamma function, defined as "^{z) = Y{z)^ where z is 
a complex number with a positive real component, F is 
the gamma function, F' its derivative, and ac = ^r] = 
§^ mbfib-mafia ^ Equation ^ gives Tc as a function of 
the mass and chemical potential asymmetry, encoded in 
77. Due to the highly non-linear term 3^{ac 7^ 0) in this 
equation, it is possible to find analytical solutions for 
Tc only when the Fermi momentum of the two particle 
species match, Pp = Pp, where = ^/2rn^/I^. This 
can happen when = 7715 and /i^ = which results 
in 3^{ac = 0) = -2j - 41n(2), and a = 1/2, yielding 
the standard BCS result Tc/Aq = ^ or, in a very exotic 
situation, where 7^ and jj^a 7^ /i6 but, in spite 
of the asymmetries, the particles have the same Fermi 
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surface, m^/ia = ^6^6, giving Tc/Aq = 2a^ [H, [2^. 
As can be seen from the expression for ac, the critical 
temperature Tc given by Eq. ([5j) is a maximum when 
the Fermi surfaces for the a and b particles are equal, 
since 3^(ac = 0) is a minimum. We obtain numerically 
the critical temperature of an imbalanced system, not 
restricted to the previously analyzed equal Fermi surfaces 
cases. First we note that Eq. ([5]) can be written in terms 
of dimensionless quantities as 
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\27rTc/AoJ 2 V 



1 /Sm /i S^\\ 



where Sm = rrih — rua^ rh = mb -\- ma, /i = /ia + 
Sjii = jih — l-ia- Since we also consider mass asym- 
metry, it is convenient to rewrite <j and M in terms 
of the nondimensional mass ratio — 1 < ^ < 1, as 



^ = ^ ^^'^M = \ma{l+'-^). As the 

chemical potential difference dji increases, Eq. ([6]) has 
two solutions, i.e., there are two temperatures at which 
A = 0. This can be seen in Fig. ([1]), where we plot 
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for different values of J/i/Ao, with equal masses, as a 
function of T. Throughout the paper we use, when nec- 
essary, /i/Ao = 7. We can see that, for 6fi/Ao small 
enough there is only one solution of Eq. (|6]). At a certain 
value of the chemical potential asymmetry J/io, Eq. (|6|) 
starts to develop two solutions, Tc,i and Tc,2. We find, 
as we shall see below, that Tc,2 corresponds to a second 
order phase transition temperature, whereas Tc,i corre- 
sponds to the temperature at which the local maximum 
of the thermodynamic potential disappears. Increasing 
the chemical potential asymmetry beyond Sfio a criti- 
cal value is reached Sfic^ after which the system presents 
a first order phase transition at a critical temperature 
Tc. Finally, at a higher chemical potential asymmetry 
^l^cc/ Aq = >/2, the Chandrasekhar-Clogston limit of 
fermionic superfluidity [27, 28], the system is a normal 
Fermi liquid, at any temperature, even T = 0. 

The critical temperature as a function of Sjj. is pre- 
sented in Fig. ([2|) for equal masses. Curves I and III 
are the second and first order phase transition lines, 
respectively. Curves II and IV are the unstable and 
met ast able, second and first order, respectively, phase 
transition lines. Below the tricritical point one has to 
properly minimize the free energy rather than using the 
gap equation as the transition becomes first order. Thus, 
the second order transition line comes simply from the 
solution of the gap equation. To construct the first or- 
der line we numerically search for the temperature Tc at 
which Q{Amin^Tc) = Q{A = 0,Tc), where Amin is the 
non-trivial minimum of Q. 

Along curves II and IV A = 0, but does not corre- 
spond to the absolute minimum of Q. Thus the true sta- 




FIG. 2: The critical temperature as a function of the chem- 
ical potential asymmetry, plotted from Eq. (|5]), in the weak 
coupling approximation. 



ble phases of an imbalanced fermionic system are below 
curves I and III (the superfluid phase) and above them 
(the normal phase). The tricritical point (^/ic,^c) = 
(I.23O8A0, 0.2713A0), and the chemical potential differ- 
ences Sjj.0 and Sjucc are shown in the figure. The small 
bumps at line IV are due to a not high enough numerical 
accuracy, but do not affect the qualitative behavior of 
this curve. In Fig. (|3]) we show the critical temperature 
as a function of Sjii for different (positive) mass ratios 
{Sm/m runs from (outer curve) to 0.25 (inner curve), 
with fi /Ao 2). 

We can derive analytical expressions for d/j^o and S/j^c 
even when Sm 7^ 0. Looking at Fig. ([1]) we see that 
the chemical potential asymmetry at which Eq. (|6]) be- 
gins to present two solutions, 6fio^ occurs when f{T = 
0) = 0. If we observe that, as T ^ 0, ^{ac) goes ex- 
actly to 2/n(|ac|/7r), then Eq. (j6|) can be rewritten as 

in{^T^) = 52^n(|is^(^^ + |^)|). This 
implies in 



Ao 



±2a- 



Sm ji 
m Ao 



(8) 



For 5m = (a = 1/2) the two critical temperatures start 
to show up at = ±1, exactly as obtained numeri- 
cally in Fig. ([2j) for > 0. Conversely, for fixed chem- 
ical potential imbalance and (artificially) varying mass 
asymmetry, the two Tc will begin to appear at 
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The behavior of Tc/ Aq as a function of ^ is similar 
to that presented in Figs. ([2]) and ([9|). The tricritical 
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point is obtained by imposing the vanishing of 

the quadratic (= A) and the quartic (= B) terms in the 
Landau expansion of the free energy, Eq. ([2j). The co- 
efficients A and B are given by f(T) and i<idf{T)/dT^ 
respectively, where is a constant [29|. Graphicahy this 
can be seen in Fig. ([T]) for Sfi = Sjic^ with Sfic correspond- 
ing to the chemical potential asymmetry where Tc occurs 
at the minimum of f{T) and above which there are no 
more solutions for /(T) = 0. From Eq. ([7j) we obtain 



df{T) 
dT 



1 

T 



2T 



We can rewrite 



S"(a) = 0. 



(10) 



as a 



2 rua+mb 

^(^A^ + A^)- The solutions of Eq. i^,aT{a) = 2, 
are ±ao, where ao = 1.1343, i.e., ao is the argument that 
minimizes /(T). Thus, the temperature that minimizes 

fiT) is T,,^,, = A^(|^^ + ^|).As discussed above, 
the (tri-)critical chemical potential Sfic is determined by 

Tc,min Tc m Eq. (O, ^(^^^^ + ^) ^c,mm 

^g-i9^(ao)^ resulting in 



(11) 



where a = cr(^). From the equation above one clearly 
sees that a positive mass imbalance reduces the critical 
chemical potential imbalance. Computing ^ for Sm = 
(g-i9^(ao) ^ 3.4089), we find ^ ^ ±1.2308. Also, when 

Sm = 0, we can use Eq. ([5]) and replace ^e~^^^^°^ = 
2Tc/Ao in Eq.([TI]), obtaining 



Ao 



= ±4ao-^, 

^0 



(12) 



and the critical temperature for the equal masses system 
at the critical chemical potential results ~ 0.2713, 
which compares exactly with those numerically obtained 
in Fig. ([2]). This result should be contrasted to that of 
the standard (symmetric) BCS result = ~ ^ 0.5669. 
The point {Sfic^Tc) is a tricritical point, with a line of 
first-order transition emerging, and hitting the Sjii axis 
at (5/icc,0). Note that the phase diagram of the equal 
mass system is symmetric with respect to the Tc /Ao 
axis. From Eq. (|lip . one can see that for ^ > ^ 



and A 



^6m ^ where ^ = , 

^e~2^(^o) ~ 1.2312, the phase diagram can be en- 
tirely shifted to the left (right) of the Tc/ Ao axis, ex- 
hibiting two negative (positive) values for the tricritical 
chemical potential imbalance, which we denote Sfici/^o^ 
and 6fic2/^0' Beyond these tricritical points there are 
also generally the CC points limiting superfiuidity be- 
tween ((^/icci/Ao, 0) and ((5/ic'C2/Ao, 0). As we will dis- 
cuss later, in trapped systems the chemical potentials 
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FIG. 3: The critical temperature as a function of S/j^/Aq for 
several mass ratios 6m /m, as indicated in the figure, with 
m/Ao = 2. 



of the two species vary with the position from the cen- 
ter to the edge of the trap. This can assure the varia- 
tion of Sfi/Aci necessary to the formation of shell struc- 
tures [3, [ill, [lO, IIH, composed of a superfiuid phase 
{Sficci/^o < ^m/Ao < Sficc2/^o) sandwiched between 
a normal core {6fi/Ao < Sficci/^o) and an outer normal 
phase (^/i/Ao > 811002/ ^q)^ in real space. Thus we see 
that only mass imbalanced trapped systems allow shell 
structure configuration. Physically, one can understand 
the displacement of the (weak coupling) phase diagram 
analysing the position of the top of the dome: the maxi- 
mum critical temperature corresponds, for a given mass 
asymmetry, to the chemical potential difference at which 
the Fermi surfaces match. It is an easy task to show that 
the maximum Tc in Eq. ([5]) occurs when ^ = — or, 
equivalently, m^/ia = ^6/^6 (see also Eq. (j6j)). Note from 
this expression that a higher positive mass asymmetry 
implies in a more negative chemical potential asymmetry 
at the top of the dome. 

We can proceed in the same way as in the case of vary- 
ing and find the (tri-)critical mass imbalance 
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with Sfi fixed. Equivalently, when S/i/Aq = the phase 
diagram Tc /Ao vs. 5m /in is symmetric with respect to 
the 5m /in = axis, and the dome is shifted to the left 
or to the right depending on the sign of a nonvanishing 
^/i/Ao. 

Numerical Solutions, a. Weak coupling. To find the 
gap which correspond to stable solutions of the system, 
as well as the first order, second order and the metastable 
phase transitions, we numerically implement the thermo- 
dynamic potential expression, Eq. ([2]), and its derivative 
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FIG. 4: The plot oi G = as a function of A/Ao for 

various temperatures. The zeros of G give the gaps that are 
candidates for the minimum of Q. For T greater than ^ .41 Aq 
there are no more solutions for G, which means that this is 
the actual critical temperature of the system (see Fig. (|2])). 





6n=1.lA„ 
X k^a= -0.5 
\ T=0.06a^ 

\ 


/ 

/ 
/ 

/ 

/ 

/ 




\ 


/ 









FIG. 5: Q/At as a function of A/A^ for kpa = -0.5, 
S/j. = l.lAo, and T = O.OGAq. For these parameters the gap 
equation has two solutions, as shown in Fig. (|4]): the smaller 
gap is related to a local maximum of Q, and the larger gap 
corresponds to a minimum of Q. 



with respect to A, and solve them for several masses 
and chemical potential asymmetries, at different tem- 
peratures. The derivative of Q is proportional to 2A, 
so A = is always a extremum of Q. To focus on 
the A 7^ solution of dQ/dA = 0, we introduce the 
(nondimensional) function G = ^^Z"^* Thus, the ze- 
ros of G give the non-trivial extrema of Q. We plot, in 
Fig. (jl]), G as a function of A/Aq. The coupling g can 
be related to the two-body s-wave scattering length a by 
- = -o^ + r 7#T3 where 6c = t^- Then we use 
this expression to trade g for the dimensionless parameter 
kpci and, doing that, we can establish contact with ex- 
periment. The zeros of G for a given set of parameters /i. 



Sfi, M, T, and kpa give the gaps that are candidates for 
the minimum of Q. To avoid the consideration of unsta- 
ble solutions [l^ , we also verify the behavior of Q versus 
A with these same parameters. We take a > (5/io 
and a T < Tc,i since the system would be in the normal 
phase with these parameters (see Fig. ([2]), that has been 
plotted from Eq. ([5])). We can observe that there are 
two gaps for Sfi = I.IAq, and T = O.OGAq. The smaller 
gap (= Ai) does not correspond to a minimum of ft, it 
is a local maximum, whereas the larger gap (= A2) is 
related to the (superfluid) stable solution, as can be seen 
in Fig. ([5]) where we plot as a function of A/Aq. 

Still in Fig. (|4j), we see that as the temperature is in- 
creased from T = 0.06 Ao up to T = 0.12Ao which is 
~ Tc,i the system ceases to have two Tc. This is an indi- 
cation that the lower Tc of Fig. ([2]) (remember, the two 
Tc appear only for Sfio < 6fi < Sfic) is the critical tem- 
perature of the unstable gap. Another way of seeing that 
the smaller gap corresponds to an unstable phase is the 
observation of the quasiparticle excitations depicted in 
Fig. ([8]). With the smaller gap the quasiparticles of the 
lower branch would have negative excitations, exhibit- 
ing the same unstable zero temperature behavior 0, 0] • 
For Tc,i < T < Tc,2 there is only one solution, which 
always corresponds to a stable phase. For even greater 
temperatures, or T > Tc,2 there are no more solutions for 
the (second order transition) gap equation, Eq. (|4]). For 
Sjii > djic we compare the thermodynamical potential 
evaluated at A = with the same quantity evaluated 
at the solutions Ao(T) of the equation G = 0, which 
gives the extrema of the thermodynamic potential. The 
temperature where Vt{A = 0,T) = ^](Ao(T),T) corre- 
sponds to the first order phase transition critical temper- 
ature Tc- For temperatures greater than Tc, Eq. G = 
can still present solutions, but they do not correspond to 
the absolute minimum of Vt. Curves IV (associated with 
metastable phases) in figures ([2]) and ([9]) correspond to 
the temperatures above which G = has no more solu- 
tions. 

Fig. shows the results for the dependence of the 
critical temperatures with the mass asymmetry. As an- 
ticipated, the behavior is quite similar to that presented 
for the chemical potential asymmetric case, as discussed 
above. Nevertheless, we must remark that the depen- 
dence of the critical temperature Tc with mass asymme- 
try is not the same as for the chemical potential asymmet- 
ric case, as can be viewed from Eq. ([6]), since <j depends 
on ^ but not on dji. In obtaining this figure, we kept 
rua constant (= Aq) and varied the mass of particle h. In 
experiments, it would be equivalent to measure the crit- 
ical temperatures of the system with species a fixed, for 
several h particles with different masses, but always keep- 
ing the chemical potential asymmetry fixed. Above some 
specific mass asymmetry, Smo/m^ the system presents 
two critical temperatures (as happens for varying 
with fixed Sm), one corresponding to the temperature 
where a local maximum of disappears (line II), and 
the other to the critical temperature where the system 
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FIG. 6: The critical temperature as a function of the (posi- 
tive) masses asymmetry for ^/i = 0.4Ao, in the weak couphng 
hmit. 



undergoes a second order phase transition (Hne I). The 
behavior of the system as a function of the mass asym- 
metry also possesses a tricritical point Tcp, where the 
second order phase transition line (I) meets the first or- 
der phase transition line (III) at 5mc/rh^ as we already 
discussed. There is also a critical mass asymmetry, above 
which the system no longer presents a superfluid phase, 
even at T = 0. This would be the equivalent of the 
Chandrasekhar-Clogston (CC) limit of superfluidity for 
the case of mass imbalance. Besides, for temperatures 
below those of curve IV, the system could still present a 
(met ast able) local minimum of Vt. 

It is interesting to note that the possibility of pairing 
of particles with different masses and chemical potentials 
results in situations not present in the purely chemical 
potential {5m = 0) or purely mass {S/j. = 0) asymmetric 
cases. If we have, for example, > m5 (which implies 
in a negative ^) in considering /i^ > /i^, the phase di- 
agram behaves as depicted in Fig ([7j), i.e., the critical 
temperature is higher for a non vanishing mass asymme- 
try. We can observe it, for instance, in the analytical 
solution for the weak coupling regime, as a consequence 
of the fact that Tc given by Eq. (|6|) does not have its max- 
imum value for 5m /in = 0, but rather when the Fermi 
surfaces match. 

b. Moderate to strong coupling. Now we investigate 
the solutions of Eq. @ in the moderate to strong cou- 
pling limit, kpa = —10, and show the resulting Tc/Aq 
as a function of J/i/Ao in Fig. (|9|). As one can sees, this 
curve has the same qualitative behavior of that in the 
weak coupling limit. We obtain next, the solutions of the 
gap equation, the thermodynamic potential and the quasi 
particle excitations, which are shown in Figs. (p!Q|) to (p!2|) 
below, respectively. As seen in the weak-coupling regime 
(Fig. ([8])), for some temperatures (T < Tc,i, where Tc,i 
is the temperature at which the unstable gap vanishes, 
for any kpci) the quasiparticles of the lower branch have 



FIG. 7: The critical temperature as a function of the (nega- 
tive) masses asymmetry for 6/^ = 0.4Ao, in the weak coupling 
limit. The critical temperature presents its maximum value 
at a non vanishing mass asymmetry. This phase diagram is 
the continuation of that of Fig. (|6|) for the negative values of 
Sm/rfi. The higher Tc is located in the top of the dome when 
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FIG. 8: The quasiparticles excitation energies behavior as a 
function of the two gaps Ai < A2 found for T = O.O6A0. 
With the smaller gap, the lower branch has quasiparticles 
{^k,b) with negative energy. 



negative excitations (Fig. (fT2|)) between two momentum 
values, and consequently, two Fermi surfaces with gap- 
less modes (BP2). This is in the heart of the instability 
of the breached-pair or phase separation in momentum 
space phases both at zero [6| and finite temperature, as 
we have shown here. However, it is worth to mention that 
a state where the lower branch excitation crosses zero 
once, defi ning on ly o ne Fermi surface (BPl), is found to 
be stable ^l[Ii,[32[. 

To take into account the trap effects, usually one em- 
ploys the local density approximation (LDA) via jiia 
^(Oa = Ma ~ "^(^5 where /i^ are the (global) chemical 
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FIG. 9: Tc/ Ao as a function of S/j./ Aq for Sm — from a 
numerical evaluation of Eq. (|4]). This curve has the same 
qualitative behavior of that shown in Fig. (|2]), obtained in 
the weak coupling limit, but clearly shows that diicikpa — 
-10) > 6^c(kFa = -0.5) . 



FIG. 11: VL/At as a function of A/A^ for kpa = -10, 
5 11 = l.lAo, and T = O.OGAq. As observed in the weak 
coupling regime (see Fig. (|5])), for these parameters the gap 
equation has two solutions, as shown in Fig. (|10|) : the smaller 
gap is related to a local maximum of and the larger gap 
corresponds to a minimum of Q. 
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FIG. 10: The plot of G = ^3^^ as a function of A/Ao for 

various temperatures. The zeros of G give the gaps that are 
candidates for the minimum of Vi. The Tc in the strong cou- 
pling limit, ~ .44Ao is larger than that of the weak-coupling 
limit, for same ^/i, (see Fig. (|4|). 



potentials introduced before, and V(f) is the trapping 
potential. As we mentioned already, trapped fermions 
with unequal masses can exhibit a shell structure. It is 
very interesting to note that our MF weak-coupling ana- 
lytical results for = in Eq. (|6]), are independent of 
/i(r) = /i — 2V(r), and depends only on ^/i(r) = Since 
the shell structure manifests only for trapped fermions 

with unequal masses [3, m, Hi, this is an indication 

that the dependence of the main structure of the zero and 
finite temperature phase diagram on /i(r) is main- 
tained in the strong coupling limit. 

Discussion and Conclusion. We have investigated 
temperature effects and thermodynamical phase transi- 



tions in fermionic gases composed by two particle species 
whose Fermi surfaces or densities do not match. We have 
observed two of the three phase transitions we mentioned 
in the introduction. At zero temperature (and at low 
fixed T, as in Figs. ([5]) and (fTT]) ), increasing the chemical 
potential asymmetry, at Sjiic the minimum of ^^/Aq jumps 
from A/Ao = 1 (A/Aq = 1.1 in Fig. ([II])) to A/Aq = 0. 
At fixed d/ii < dfic-, increasing T the minimum of ^^/Aq 
goes smoothly from some A/Aq < 1 to A/Aq = 0. For 
5 lie < dji < djicCi increasing T the minimum of ^^/Aq 
jumps from some A/Aq < 1 to A/Aq = 0. We observed 
the same conclusions for the varying mass situation. Be- 
sides the second and first order lines, respectively, asso- 
ciated with the phase transitions we just mentioned, we 
have found the unstable and metastable lines related with 
these phase transitions. We have shown that for specific 
values of the chemical potential (and/or mass) asymme- 
try the thermal gap equation of imbalanced systems has 
two solutions. The smaller gap, corresponding to the BP 
phase, always represents unstable solutions. We have 
found for the first time, in the weak coupling limit, an 
analytical expression for the tricritical chemical potential 
difference, which in turn depends on the mass asymme- 
try. A positive (negative) mass asymmetry reduces (in- 
creases) the positive tricritical chemical potential imbal- 
ance, as a consequence of the displacement of the dome to 
the negative (positive) side of the Tc/ Aq axis. These con- 
clusions were obtained also in the strong coupling regime. 
We have also shown that raising the interaction param- 
eter kpa-, results in a small increasing of the (critical) 
chemical potential differences Sfic and Sficc and, equiva- 
lent ly, of the critical mass differences. We also have iden- 
tified the presence of metastable states beyond the first- 
order transition lines (curves III in Figs. ([2|), ([6|), ([7j), and 
([9|)) and the temperatures at which these metastable so- 
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lutions disappear in a first order phase transition (curves 
IV in the same figures). Given the current possibility 
of controlling several parameters in ultracold fermionic 
systems, such as the interaction strength, densities, trap- 
ping potential, and temperature, it would be possible the 
experimental preparation of long-lived metastable states 
of imbalanced superfiuids, allowing the study of thermo- 
dynamic and decay process of these states, a subject of 
interest in several areas of physics [sl, HI, [35[ . We hope 
that our results for the imbalanced masses can be veri- 
fied, considering the limitations of the MF theory, in the 
new generation of experiments that could in principle be 
set up with two kinds of fermionic atoms with opposite 
spins and different masses (36| . 



FIG. 12: The quasiparticles excitation energies behavior as a 
function of the two gaps Ai < A2 found for for kpa = —10, 
and T — O.O6A0. As in the weak coupling limit (Fig. 6), with 
the smaller gap, the lower branch has quasiparticles (Sfc,fe) 
with negative energy. However, in the strong coupling regime 
the range of momenta at which the particles have negative 
energy (meaning that they are single) is bigger than in the 
weak coupling limit. 
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